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ABSTRACT 

Hughston has recently proposed a stochastic extension of the Schrodinger 
equation, expressed as a stochastic differential equation on projective Hilbert 
space. We derive new projective Hilbert space identities, which we use to 
give a general proof that Hughston's equation leads to state vector collapse 
to energy eigenstates, with collapse probabilities given by the quantum me- 
chanical probabihties computed from the initial state. We discuss the relation 
of Hughston's equation to earlier work on norm-preserving stochastic equa- 
tions, and show that Hughston's equation can be written as a manifestly 
unitary stochastic evolution equation for the pure state density matrix. We 
discuss the behavior of systems constructed as direct products of indepen- 
dent subsystems, and briefly address the question of whether an energy-based 
approach, such as Hughston's, suffices to give an objective interpretation of 
the measurement process in quantum mechanics. 
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I. INTRODUCTION 

A substantial body of work [1] has addressed the problem of state vector collapse 
by proposing that the Schrodinger equation be modified to include a stochastic process, 
presumably arising from physics at a deeper level, that drives the collapse process. In 
particular, Gisin [2], Percival [3], and Ghirardi, Pearle, and Rimini [4] have constructed 
equations that preserve the norm of the state vector, which in the approximation that the 
usual Schrodinger Hamiltonian dynamics is neglected are shown [4] to lead to state vector 
collapse with the correct quantum mechanical probabilities. An alternative approach to con- 
structing a stochastic extension of the Schrodinger equation has been pursued by Hughston 
[5] , based on the proposal of a number of authors [6] to rewrite the Schrodinger equation as 
an equivalent dynamics on projective Hilbert space, i.e., on the space of rays, a formulation 
in which the imposition of a state vector normalization condition is not needed. Within 
this framework, Hughston [5] has proposed a simple stochastic extension of the Schrodinger 
equation, constructed solely from the Hamiltonian function, and has shown that his equa- 
tion leads to state vector reduction to an energy eigenstate, with energy conservation in the 
mean throughout the reduction process. In the simplest spin-1/2 case, Hughston exhibits 
an explicit solution that shows that his equation leads to collapse with the correct quan- 
tum mechanical probabilities, but the issue of collapse probabilities in the general case has 
remained open. 

Our purpose in this paper is to further investigate the structure and properties of 
Hughston's equation, proceeding from new identities in projective Hilbert space derived in 
Sec. II. A principal result will be the proof in Sec. Ill (using the martingale or "gambler's 
ruin" argument pioneered by Pearle [7]) that in the generic case, with no approximations, 
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Hughston's equation leads to state vector collapse to energy eigenstates with the correct 
quantum mechanical probabilities. The relation of Hughston's equation to earlier work on 
norm-preserving equations is discussed in Sec. IV, and the density matrix form of Hughston's 
equation, which gives a manifestly unitary stochastic evolution on pure states, is given in 
Sec. V. In Sec. VI we examine the stochastic evolution of an initial state that is constructed 
as the product of independent subsystem states. Finally, in Sec. VII we discuss whether 
an energy-based approach to stochastic evolution (as opposed to approaches [8] based on 
spontaneous localization) suffices to give a satisfactory objective description of the evolution 
of a state during the quantum mechanical measurement process. 

II. PROJECTIVE HILBERT SPACE AND SOME IDENTITIES 

We begin by explaining the basic elements of projective Hilbert space needed to un- 
derstand Hughston's equation, working in an n -|- 1 dimensional Hilbert space. We denote 
the general state vector in this space by l^;), with z a shorthand for the complex projections 
2;°, z^, of the state vector on an arbitrary fixed basis. Letting F be an arbitrary Hermi- 
tian operator, and using the summation convention that repeated indices are summed over 
their range, we define 



so that (F) is the expectation of the operator F in the state l^;), independent of the ray 
representative and normalization chosen for this state. Note that in this notation (F^) and 
{FY are not the same; their difference is in fact the variance [AF]^, 




(la) 



[AF]' = (F") - {Ff 



(16) 
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We shall use two other parameterizations for the state \z) in what follows. Since (F) is 
homogeneous of degree zero in both 2;" and let us define new complex coordinates P by 

p=z^/^^ ?=^V^o, j = l,...,n, (2) 

which are well-defined over all states for which z'^ ^ Q [9]. Next, it is convenient to split 
each of the complex numbers P into its real and imaginary part t-Jj, and to introduce 
a 2n component real vector x", a — 1, ...,2n defined by — t^, — t\, — t\, x'^ — 
t^, ...,x^"'~^ = t^, x'^"' = P}. Clearly, specifying the projective coordinates P or x"" uniquely 
determines the unit ray containing the unnormalized state | z) , while leaving the normaliza- 
tion and ray representative of the state l^;) unspecified. 

As discussed in Refs. [6], projective Hilbert space is also a Riemannian space with 
respect to the Fubini- Study metric Qapi defined by the line element 

\{z\z + dz)\'^ 



{z\z) {z + dz\z + dz) 
Abbreviating z'^ = z-z, a simple calculation gives 



(3a) 



d d 

gap = 4(5a^z • z - z''z^)/{z ■ zf = 4— — logz • z . (36) 

Because of the homogeneity conditions z'^gap — z^gajS — 0, the metric gap is not invertible, 
but if we hold the coordinates z^ fixed in the variation contained in Eq. (3a) and go over 
to the projective coordinates P , we can rewrite the line element of Eq. (3a) as 

ds^ = gjkdt'dt'^ , (4a) 

with the invertible metric [9] 

m+t't%,-pt] 

93k- (l+ri-)2 ' ^ ' 



with inverse 

gJk^^^l+m^Sjk + tn") . (4c) 

Reexpressing the complex projective coordinates P in terms of the real coordinates x"", the 

line element can be rewritten as 

ds^ —Qabdx'^dx^ , 

_ 4[(1 + a; - {x^x^ + uJacX^oObdX'^)] 

g""^ =^{1 + X^X^){Sab + X^X^ + UJacX'^UJbdX'^) . 

Here ujab is a numerical tensor whose only nonvanishing elements are 

uja=2j-i b=2j = 1 and u!a=2j b=2j-\ = —1 for j = 1, ...,n. As discussed by Hughston, 
one can define a complex structure Jj' over the entire projective Hilbert space for which 
JaJb'^Ocd = gab, JaJb" = -^a> ^uch that Viab = QbcJa ^nd Qf'^ = g'^'^jj' are antisymmetric 



tensors. At x = 0, the metric and complex structure take the values 
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gab =45a6 , g"^ = -5ab 



Ja —^ab , ^ab — ^^ab , — '^^ab 



(5) 



Returning to Eq. (la), we shall now derive some identities that are central to what 
follows. Differentiating Eq. (la) with respect to z", with respect to z^, and with respect to 
both z°' and z^, we get 

{^\^)^ =Fa,z^ - iF)z'^ ' 

{z\z)^^z^F^f,-{F)z^ , (6a) 

HAz)[Fa/3 - 5a(,(F)] + 2z^z^{F) - z^F.^z^ - z^F^^z^ . 
Writing similar expressions for a second operator expectation (G), contracting in various 

combinations with the relations of Eq. (6a), and using the homogeneity conditions 

dfz'' dzp dz'^dzp dz^'dzP ^ ' 
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to eliminate derivatives with respect to 2;°, we get the following identities, 

(FGF) - {F^){G) - {F){FG+GF) + 2(f )=(G) = =2V(F)V'(F)V,V4G), 

(7a) 

with Va the covariant derivative constructed using the Fubini-Study metric affine connection. 
It is not necessary to use the detailed form of this affine connection to verify the right 
hand equalities in these identities, because since (G) is a Riemannian scalar, VaVfe(G) = 
Vadb{G), and since projective Hilbert space is a homogeneous manifold, it suffices to verify 
the identities at the single point x = 0, where the affine connection vanishes and thus 
VaVb(G') = dadb{G). Using Eqs. (7a) and the chain rule we also find 

-Va[{F^)-{Ff]V%G) = -l-iF^G + GF^) + {F^){G) + {F){FG + GF)-2{Fy{G) . (76) 

When combined with the final identity in Eq. (7a) this gives 

D =V%F)V\F)VaV,{G) - ^Wa[{F') - {Ff]W^{G) 

=^{2FGF - F^G - GF'^) (7c) 
= -l([F,[F,G]]) , 

with [ , ] denoting the commutator, from which we see that D vanishes when the operators 
F and G commute. 

An alternative derivation of Eq. (7c) proceeds from the fact, noted by Hughston, 
that (for F self-adjoint) 

= Q«^V6(F) (8a) 



is a Killing vector obeying 

VeCF + V^ecF^O . (86) 

Using the identity (F^) - {Ff = Vb{F)V^{F), which is the F = G case of the middle 
equality of Eq. (7a), we rewrite D of Eq. (7c) as 

D = V''{F)V\F)VaVb{G) - V\F)V%G)VaVb{F) . (9a) 

This can be rewritten, using the identity fl'^'Tlcb — (^c> antisymmetry of fl, the fact that 
Q commutes with the covariant derivatives, and the Killing vector definition of Eq. (8a), as 

(96) 

We now use the Killing vector identity of Eq. (8b) on the final factor in each term, giving 

D = CcF^eFn'^'S/'U - ^cF^eG^'^VXaF ■ (9c) 

Exchanging the labels e and c in the first term, and exchanging the labels a and e in the 
second term, we get 



(10a) 



D HcFieF^^^VXaG + icFiaG^^'V XeF 
HcFV^^i^^'^eFU] ■ 

Substituting the Killing vector definition of Eq. (8a), this becomes 

D =Q,hV\F)V'[Q'''QefV^{F)QagV%G)] 

^ncbV\F)V'[ngfV'{G)Vf{F)] (106) 
= -i([F,[F,G]]) , 

where to get the final line we have twice used the first identity in Eq. (7a). This completes 
our geometric derivation of Eq. (7c). 



III. HUGHSTON'S EQUATION AND STATE VECTOR 
COLLAPSE PROBABILITIES 



Let us now turn to Hughston's stochastic differential equation, which reads 

dx" = [2Q"''V6(//) - ^a'^W''V]dt + (7V"(i/)dWt , (11a) 

with Wf a Brownian motion or Wiener process, with a a parameter governing the strength 
of the stochastic terms, with H the Hamiltonian operator and {H) its expectation, and with 
V the variance of the Hamiltonian, 

V = [AH]^ = (H^) - {Hf . (116) 

When the parameter a is zero, Eq. (11a) is just [6] the transcription of the Schrodinger 
equation to projective Hilbert space. For the time evolution of a general function G[x\, we 
get by Taylor expanding G[x + dx] and using the Ito stochastic calculus rules [10] 

[dWtf = dt , [dtf = dtdWt = , (12a) 

the corresponding stochastic differential equation 

dG[x\ = ndt + aV aG[x\V\H)dWt , (126) 

with the drift term ^ given by 

11 = 2^l'''VaG[x\Vb{H) - ^^a^V^VVaGix] + ]^a''V''{H)V\H)V a^bGix] . (12c) 

Hughston shows that with the cr^ part of the drift term chosen as in Eq. (11a), the drift term 
11 in Eq. (12b) vanishes for the special case G[x\ = (if), guaranteeing conservation of the 
expectation of the energy with respect to the stochastic evolution of Eq. (11a). But referring 
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to Eq. (7c) and the first identity in Eq. (7a), we see that in fact a much stronger result is also 
true, namely that /i vanishes [and thus the stochastic process of Eq. (12b) is a martingale] 
whenever G[x] — (G), with G any operator that commutes with the Hamiltonian H. 

Let us now make two applications of this fact. First, taking G[x] = V = {H"^) — {H)"^, 
we see that the contribution from (H^) to // vanishes, so the drift term comes entirely from 
— {Hy. Substituted this into /i gives —2{H) times the drift term produced by {H), which is 
again zero, plus an extra term 

-a^V''{H)V\H)Va{H)Vb{H) = -aV^ , (13a) 

where we have used the relation V — Va{H)V"'{H) which follows from the F — G — H 
case of the middle identity of Eq. (7a). Thus the variance V of the Hamiltonian satisfies the 
stochastic differential equation, derived by Hughston by a more complicated method, 

dV = -a'^V^dt + aVaVV''(H)dWt . (136) 

This imphes that the expectation E[V] with respect to the stochastic process obeys 

E[V{t)] = E[V{0)] - f dsE[V{sf] , (13c) 

Jo 

which using the inequahty < E[{V — E\y]Y] — E\y'^] — E\yY gives the inequahty 

E[V{t)\ < E[V{0)] - f dsE[V{s)Y . (13d) 

Jo 

Since V is necessarily positive, Eq. (13d) implies that i?[y(cxD)] = 0, and again using positiv- 
ity of V this implies that V{s) vanishes as s — > oo, apart from a set of outcomes of probability 
measure zero. Thus, as concluded by Hughston, the stochastic term in his equation drives 
the system, as f — > oo, to an energy eigenstate. 
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As our second application of the vanishing of the drift term ji for expectations of 
operators that commute with i7, let us consider the projectors He = |e)(e| on a complete 
set of energy eigenstates |e). By definition, these projectors all commute with H, and so the 
drift term fi vanishes in the stochastic differential equation for G[x\ = (He), and consequently 
the expectations £'[(ne)] are time independent; additionally, by completeness of the states 
|e), we have ^^i^e) — 1- But these are just the conditions for Pearle's [7] gambler's ruin 
argument to apply. At time zero, i5[(ne)] = (He) = Pe is the absolute value squared of the 
quantum mechanical amplitude to find the initial state in energy eigenstate |e). At i = oo, 
the system always evolves to an energy eigenstate, with the eigenstate |/) occurring with 
some probability Pf. The expectation £'[(ne)], evaluated at infinite time, is then 

E[{Ue)]=lxPe + J2^xPf = Pe ] (14) 

hence Pe — Pe for each e and the state collapses into energy eigenstates at t — oo with 
probabilities given by the usual quantum mechanical rule applied to the initial wave function. 
This conclusion clearly generahzes to the stochastic equation 

j 3 

with the Hj a set of mutually commuting self- adjoint operators that commute with if, with 
Vj = (Hj) — {HjY, and with the dWl independent Wiener processes obeying dW^dW^ — 
5^^dt. Following the same method used in obtaining Eq. (13b), and defining C^j = V a{Hk)V{Hj) = 
{HkHj)-{Hk){Hj), one finds 

(^14 = -(T'^C2^.dt + aVaVfc^V"(if,)rfM// , (156) 

j 3 
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and therefore 

E[Vk{t)] = EiVkiO)] - a' f dsTEiCkjisf] . (15c) 

Jo . 

Since E[Clj] > E[Ckj]'^, we have as before 

E[Vk{t)] < E[Vkm - r ds V E[Ckj(s)]' , (15d) 

Jo ^ 

which imphes that each E[Ckj{s)] approaches zero as s — > oo. Hence for each k,j we have 
at large times 

E[{H,H,)] - E[{H,){H,)] ^ , (15e) 

and so there is an effective quantum decorrelation of commuting observables. Moreover, for 
k — j Eq. (15e) imphes that at large times E[Vk\ 0, which since T4 is nonnegative implies 
that Vk approaches zero apart from a set of outcomes of probability measure zero, and so 
the state evolves to a simultaneous eigenstate of all the commuting observables entering the 
process of Eq. (15a). 

IV. RELATION OF HUGHSTON'S EQUATION TO OTHER 
STOCHASTIC NORM-PRESERVING EQUATIONS 

Let us now specialize Eqs. (12b) and (12c) to the case in which G[x] is simply the 
expectation (G) of an operator G. Then by substituting Eqs. (7c) and the second equality 
in Eq. (7a), we find 

d{G) = fxdt + KdWt , (16a) 

with 

K^la[{{G,H})-2{G){H)]^^-a{{G,H-{H)}) , (166) 
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and with 

f^={-z[G,H])-^a\[H,[H,G]]) , (16c) 

where we have used { , } to denote the anticommutator. 

Let us now compare this with the evolution of (G) imphed by the stochastic state 
vector evolution 

d\z) = [adt + pdWt\\z) , (17a) 

with 

a^-iH-la'[A-{A)]' , 

(176) 

where A is a general self-adjoint operator and (^4) is defined, as in Eq. (la), by {A) = 

{z\A\z)/{z\z). For the evolution of {z\G\z), we find by the Ito rules, 
d{z\G\z) =(z|[q;^G' + Ga + (5^G(5\dt + [(3^G + G(5\dWt\z) 

^{z\ - H]dt - \a'[{G, [A - {A)f} - 2[A - {A)]G[A - {A)]]dt 

1 (18a) 

+ [A-{A)]}dWt\z) 

={z\ - i[G, H]dt - \(T^[A, [A, G]]dt + \a{G, [A - {A)]}dWt\z) . 
8 2 

When G = 1, the right hand side of Eq. (18a) vanishes, since the commutator terms vanish 
trivially and {z\A — {A)\z) = {z\A\z) — {z\z){A) = 0. Therefore the state vector evolution of 
Eqs. (17a, b) is norm preserving, and so it is consistent to choose the normalization {z\z) — 1 
in conjunction with this evolution. For general G we then have d{z\G\z) = d{G), and so 
Eq. (18a) gives an expression for d{G), which we see is identical to Eqs. (16a) - (16c) when 
the operator A is taken as the Hamiltonian H. In particular, when A = G = H we learn 
from Eq. (18a) that d{H) — aVdWt, in agreement with Eq. (12b), because 

{z\{H, H - (H)}\z) = 2[(H') - {Hf\ = 2V , (186) 
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and so the convergence argument of Eqs. (13a)-(13d) follows directly from Eq. (18a). Apart 
from minor changes in notation, the norm preserving evolution of Eqs. (17a) and (17b) is 
the one given by Gisin [2], Percival [3], and Ghirardi, Pearle, and Rimini [4], and so we see 
that this evolution is equivalent [11] to the state vector evolution in projective Hilbert space 
given by Hughston's equation. 

The evolution of Eq. (17a) can be generalized [as was done for the Hughston equation 
in Eq. (15)] to read 

d\ z) = [adt + PjdWi] \ z) , 
j 

a^-iH-la'J2[Aj-{Aj)r , (19) 
j 

with the Aj any set of mutually commuting operators. When the Aj do not all commute with 
the Hamiltonian H, it is necessary to make the approximation of neglecting the Hamiltonian 
evolution (the —iH term in a) in proving that Eq. (19) implies state vector reduction to 
the mutual eigenstates of the Aj with probabilities given by the usual quantum mechanical 
rule. Such a proof, very similar to the one given for Hughston's equation in Sec. Ill above, 
has been given by Ghirardi, Pearle, and Rimini [4]. In order to carry through the proof 
with no approximations, it is necessary to assume that the Aj are operators in the mutually 
commuting set Hj that all commute with H, as was done in Sec. III. 

V. DENSITY MATRIX EVOLUTION 

Let us now define the pure state density matrix p by 

{z\z) 
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(20a) 



in terms of which (G) is given by 



(G) = TrpG 



(206) 



Since G is a fixed operator, Eq. (18a) for d{G) can be rewritten as 
TrGdp =Trp 



i[G, H]dt - la^[A, [A, G]]dt + la{G, [A - (A)]}dWt 

o Z 



-TrG 



-i[H, p]dt - ^a'[A, [A, p]]dt + ^a{p, [A - {A)]}dWt 



(20c) 



where in the final fine we have cychcally permuted terms under the trace. Since Eq. (20c) 
holds for arbitrary self-adjoint operators G, it implies that the density matrix obeys the 
stochastic differential equation (each term of which is self-adjoint) 

dp = -i[H, p]dt - ^a'[A, [A, p]]dt + ^a{p, [A - {A)]}dWt . (21a) 

This equation can be written in an alternative form by observing that since p is a pure state 
density matrix obeying = p, we have p{A) — pTrpA — pAp. These facts imply that 

{p, [A-{A)]}=pA + Ap-2p{A) 

^p^A + Ap^-2pAp^[p,[p,A] , 

and so we can rewrite Eq. (21a) as 



(216) 



dp = -t[H, p]dt - ^a^[A, [A, p]]dt + ^a[p, [p, A]]dWt 



(21c) 



Equations (21a) and (21c) have the following properties for general A: 



(i) Since Trdp — 0, the condition Trp = 1 is preserved by the time evolution. 

(ii) After some algebra using the Ito rules, one finds that p^ — p implies that 

{p,dp} + [dp]''^dp , (21d) 
15 



which can be rewritten as [p + dpY — p + dp. Hence the evolution of Eqs. (21a, c) is 
consistent with the pure state condition. This is required by the fact that Eqs. (21a, 
c) may be derived as consequences of Eqs. (17a, b), which are a norm preserving pure 
state evolution. The condition of Eq. (2 Id) determines the coefficient of the a^...dt 
drift term in terms of the coefficient of the a...dWt stochastic term, and so the ratio 
of these two coefficients in Eq. (21c) cannot be treated as an additional adjustable 
parameter. 

(iii) Since dp — dp\ the self-adjointness of p is preserved by the time evolution. 

(iv) Time reversal invariance is violated by the stochastic terms, since when dt and i 
are reversed in sign, the term —i[H,p\dt is invariant, but the term — |(7^[yl, [vl, pjjdt 
reverses sign. 

(v) When we take the stochastic expectation of Eq. (21c), the dWt term drops out, and 

we get 



which as pointed out by Percival [3] and Ghirardi, Pearle, and Rimini [4] is a quantum 
dynamical semigroup evolution of the completely positive Lindblad [12] form. The 
stochastic expectation E[p\ is what is usually termed the density matrix; it starts 
off at i = as a pure state density matrix but then evolves, through the stochastic 
process, into a mixed state density matrix. 

(vi) The conditions for Eq. (22a) to admit stationary solutions E[p\s with dE[p\s/dt = 
are very stringent, since when the left hand side of Eq. (22a) is zero, multiplying by 



dE[p] 



i[H,E[p]]-la'[A,[A,E[p]]] 



(22a) 



dt 
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E[p\s and taking the trace gives 



0^-iTrE[p]s[H,E[p]s] 



1 



a'TrE[p]s[A,[A,E[p]s]] 



{22b) 
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Using cyclic permutation under the trace, the first term on the right hand side 
vanishes, while the second term becomes 



which can only vanish when [74,£'[p]5] = 0. Substituting this equation back into 
Eq. (22a) then further implies that [H, E[p]s] — 0. When A and H commute, 
these conditions can be satisfied with E[p]s a general function of H [see the further 
discussion of this case in (xiii) below] , but when A and H do not commute, one can 
have situations in which either E\p]s must be a multiple of the unit operator, which 
trivially commutes with both H and A, or else there are no stationary solutions and 
E[p] diverges at large times. The latter case is found in spontaneous localization 
models, as discussed for example in Section III.B.3 of Ref. [4]. 
(vii) The evolution implied by Eq. (22a) leads [13] to a monotonic increase of the von Neu- 
mann entropy, a result that can be demonstrated directly from Eq.(22a) as follows. 
Letting 



be the von Neumann (or information) entropy, we find by substituting Eq. (22a) and 
using cyclical permutation of factors under the trace, that 



(22c) 



S 



TrE[p]logE[p] 



{22d) 




dE[p\ 



[l + \ogE[p]] 



dt 



(22e) 



a'T,[\ogE\plA][A,E\p]] 
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Substituting complete sets of eigenstates \n), |m) of the nonnegative density matrix 

E[p], this becomes 
dS 1 

=g^' ^J^^log E[p]n - log E[p]m]AnmA,^n[E[p]n - E[p]m] 
n m 

n m 

(viii) Since Eqs. (21a) and (21c) are nonhnear in p, the Schrodinger dynamics described 
by them cannot be represented as an equivalent Heisenberg or dual dynamics on the 
operator G. On the other hand, Eq. (22a) is linear in p, and so as noted by Lindblad, 
the Schrodinger dynamics for E[p\ can be represented as a dual Heisenberg dynamics 
for E[G], given by 

^ = t[H, E[G]] - ^a'[A, [A, E[G]]] . {22g) 

(ix) The evolution of Eq. (21c) can be written (after some algebra, and again using 
= p) in the manifestly unitary form 

p + dp= UpU^ , U = e'^^ , (23a) 

with the infinitesimal anti-self-adjoint generator dK given by 

dK 



-iH - la^[A'^ - 2ApA, p] 



dt-^a[p,A]dWt . (236) 



Equations (21c) and (23b) thus give the stochastic unitary extension of the Lindblad 
evolution of Eq. (22a) [14]. 

Specializing to the Hughston case A — H, Eq. (21c) (which uses the pure state 
condition p^ = p) becomes 

dp = -i[H, p]dt - ^a'[H, [H, p]]dt + ^a[p, [p, H]]dWt , (24a) 
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while Eq. (22a) becomes 



dE[p] 



i[H,E[p]]-\a'[H, [H,E[p]]] 



(246) 



dt 



and the following further properties are evident: 

(x) When p = He, the projector on an energy eigenstate, then since all commutators in 
Eq. (24a) vanish we have dp — 0. 

(xi) For G commuting with H, E[d{G)] = TrCdElp] = 0, since by cyclic permutation 
inside the trace each term arising from substituting Eq. (24b) into the expectation 
of Eq. (20b) can be rearranged to have a factor [G, H]. 



(xii) For V = [AH]"^ = Ti pH^ - [Tr pH]'^, use of Eq. (24a) and the Ito calculus imply that 
E[dV/dt] = -E[R% with 



(xiii) Items (x) through (xii) are the ingredients used in Sec. Ill to prove state vector 
collapse to energy eigenstates |e) with the correct quantum mechanical probabilities 
Pe- Hence at large times, as noted by Hughston, E[p] — > '^^pelle, which explicitly 
exhibits the role of E[p] as the density matrix that evolves, under the stochastic 
process, from a pure to a mixed state form. The fact that E[p] for Hughston's 
equation approaches a stationary limit at large times is in accord with the general 
stationarity discussion given in (vi) above. 




(24c) 
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VI. BEHAVIOR OF SYSTEMS CONSTRUCTED 
FROM INDEPENDENT SUBSYSTEMS 



Let us next examine the structure of Hughston's equation for a Hilbert space con- 
structed as the direct product of independent subsystem Hilbert spaces, so that initially at 
time t — the state vector is 



We assume the Hamiltonian 



H^Y.^^ , (256) 

with acting as the unit operator on the states \zk) i k ^ i. Then a simple calculation 
shows that the expectation of the Hamiltonian (if) and its variance V are both additive over 
the subsystem Hilbert spaces, 

^ ^ ^ (25c) 

with {F()ii the expectation of the operator Fi formed according to Eq. (la) with respect to 
the subsystem wave function \z(). In addition, the Fubini-Study line element is also additive 
over the subsystem Hilbert spaces, since 

1 - ds^lA + ^ -j-r \l^Zt\zi^dZi)\^ 

{z\z){z + dz\z + dz) AJ- l^zt\zi)lyii + dzi\zi + dzt^ 

^\^\-ds\lA\^l-Yds^^l Avoids') . 
I e 

[An alternative way to see this is to use the identity log z- z — log Yle ze- z^ — log ze ■ z^ in 
Eq. (3b), along with a change of variable from z to the z^s.] As a result of Eq. (26), the metric 
g"'^ and complex structure Q"-^ block diagonalize over the independent subsystem subspaces. 
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Equations (25a)- (25c) then imply that Hughston's stochastic extension of the Schrodinger 
equation given in Eq. (11a) separates into similar equations for the subsystems, that do not 
refer to one another's x"' coordinates, but are correlated only through the common Wiener 
process dWt that appears in all of them. These correlations result in the entanglement of 
the states jz^), so that the product form of Eq. (25a) is not maintained for times t > 0, 
but subsystems \zi) already in energy eigenstates remain unentangled for all time, since the 
coefficient of dWt vanishes in their stochastic evolution equations. 

These same conclusions follow from the density matrix form of Hughston's equation 
given in Eq. (24a), in which the entanglements arising from the action of the same Wiener 
process on all subsystems are already evident, because the density matrix depends quadrat- 
ically on the normalized state vector. Considering for simplicity the case of two independent 
subsystems, substituting the t — form 

P = PiPi (27a) 

into Eq. (24a), with H = Hi + H2, we get 

dp = dpip2 + pidp2 - ^a^[Hi, pi][H2, p2]dt , (276) 

with dpi the evolution predicted by Eq. (24a) within subsystem 1, 

dpi = -i[Hi,pi]dt - la'[Hi, [Hi, p]]dt + l-a[pi, [pi, Hi]]dWt , (27c) 

o Z 

and similarly for dp2- The entangling term proportional to [Hi, pi][H2, P2]dt comes from the 
[dWt]^ contribution from the state vector evolution equation to the density matrix equation; 
it is in general nonzero, but vanishes when either [Hi, pi] = or [H2, P2] = 0, that is, when 
either of the two subsystems is in an energy eigenstate. When more than two subsystems are 
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present, the entangling term coupling pL to pe , £ Lis more complicated in structure, but 
still has a factor [Hl,Pl] and so vanishes when the subsystem L is in an energy eigenstate. 
Thus the endpoints of the stochastic evolution under Hughston's equation, which are the 
energy eigenstates, can persist indefinitely as unentangled independent subsystems in a larger 
system. 

This conclusion does not extend to the more general evolution of Eq. (21c), in which 
the stochastic process is driven by an operator A differing from the Hamiltonian, with A 
taken to be additive over subsystems. The reason is that there is now a competition between 
the stochastic terms, which are constructed from double commutators with an innermost 
commutator and the Schrodinger evolution term, which involves the commutator 

[H, p]; the stochastic terms tend to drive the system to A eigenstates, while the Schrodinger 
term coherently mixes A eigenstates, leading to evolution away from A eigenstates. Thus, 
a subsystem cannot remain indefinitely in an A eigenstate, and as a result does not persist 
indefinitely as an unentangled independent subsystem in a larger system. [These statements 
are in accord with the conclusions reached in the stationarity discussion of (vi) in Sec. V.] 

VII. DOES AN ENERGY-BASED EQUATION SUFFICE? 

In the preceding sections we have seen how Hughston's equation fits into the general 
framework of stochastic modifications of the Schrodinger equation that have been studied 
in the past. Its distinguishing feature is that the general operator A of Eqs. (17a, b) and 
(21a, c) driving the stochastic terms is chosen, in Hughston's case, to be the Hamiltonian 
H. This choice confers the advantage that the proof of reduction of the state vector to 
A eigenstates (i.e., in Hughston's case, to energy eigenstates) with the correct quantum 
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mechanical probabilities becomes exact, since it is not necessary to neglect the Hamiltonian 
evolution term. Moreover, since for Hughston's equation the stochastic expectation of the 
Hamiltonian operator E[H] is conserved in time, and since convergence to H eigenstates 
preserves the quantum mechanical predictions, any statistical test of energy conservation 
performed on the endpoint of the stochastic process will agree with the quantum mechanical 
prediction. To justify these advantages, we must now address the issue of whether an energy- 
based stochastic equation is sufficient to give an objective interpretation of state vector 
reduction [15]. 

First, we must deal with the objection that in most measurements, the quantum 
attribute being measured is not an energy; for example, in a Stern-Gerlach experiment, it 
is typically the z component of a spin. However, to perform a measurement, it is always 
necessary to couple the quantum attribute being measured to the apparatus through an 
interaction energy term Hj, in such a way that the macroscopic state of the apparatus is 
ultimately determined by the quantum attribute being measured. Thus, in the first instance, 
what is being measured is an energy, even though after amplification to macroscopic scale 
this can be converted to other forms of indication, such as pointer displacements. So from 
the point of view of the variety of quantum attributes that can be measured, Hughston's 
equation appears to be as viable as localizing approaches [8] in which A is chosen as an 
operator that produces spatial localization. 

We must next deal with the issue of whether an energy-based approach can prevent 
the occurrence of macroscopic quantum superpositions. For example, take a macroscopic 
object and displace it a macroscopic distance; the two states have the same energy, and so in 
Hughston's approach such superpositions would appear to be allowed, whereas in localizing 
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approaches they are strongly forbidden. However, this objection neglects the interactions 
of the macroscopic object with its environment, of the same type that are important in 
studies of decoherence. When such effects are taken into account, macroscopic displacement 
of a macroscopic object results in an energy shift A£^, reflecting the altered environment, 
which is sufficient, from the point of view of Hughston's equation, to lead to rapid state 
vector reduction to one displaced alternative or the other. To study this quantitatively, let 
us consider the following two environmental effects: (i) thermal energy fluctuations, and (ii) 
the surface adsorption of surrounding molecules. Hughston proposes, as have other authors 
[16], that the parameter governing the stochastic terms is of order a ~ ^pianck microscopic 
units with % — c— 1, which he shows leads to state vector reduction in a time given by 



Hence to get a reduction time of order, say, 10~^ seconds, one needs a AE' ~ 3GeV ~ 
3 nucleon masses. 

Considering first the effect of thermal fluctuations, let us consider a macroscopic 
object with N ~ 10^^ nucleon masses, so that AE ^ N2kT ^ 8GeV at room temperature 
(300° Kelvin) and AE ~ .08Gev at the 3° temperature of the cosmic microwave background. 
For such an object, thermal energy driven state vector reduction will occur in 10~^ seconds 
at room temperature and in 10~^ seconds at the temperature of the microwave background. 
Examining next the effect of adsorbed molecules, consider an object with a surface area of 
1 cm^ at room temperature in an extreme vacuum of lO^^^Torr (less [17] than the nighttime 
pressure at the surface of the moon.) Then the flux of molecules bombarding its surface is 
[17] 4 X 10® per second, so assuming a high probability for the molecules to stick, a AE of 
3GeV is attained in of order 10~® seconds, permitting a 10~® second state vector reduction 




(28) 
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time driven by the change in energy produced by surface adsorption. One can scale to other 
sizes of macroscopic object from these examples, but they suffice to show that in the normal 
range of laboratory operating conditions for measuring apparatus, environmental interactions 
produce a large enough spread of energy values to give rapid state vector reduction through 
an energy driven stochastic equation. 

^Prom a formal point of view, it is instructive to cast the above discussion of envi- 
ronmental effects in terms of the analysis of the measurement process given by Zurek [18], 
starting from Eq. (24b) for the evolution of the stochastic expectation of the density matrix. 
Zurek assumes that the total Hamiltonian H describes the system S being measured, the 
apparatus A doing the measuring, and the environment £. Thus, he writes the Hamiltonian 
as a sum of 6 terms, 

H = Hs + Hj^ + H£ + Hsa + Has + Hss , (29) 

with the first three terms giving the Hamiltonians of the system, apparatus, and environment 
in isolation from one another, and with the second three terms giving the corresponding 
interaction Hamiltonians. Zurek assumes that the interaction Hgg between system and 
environment can be neglected, and that the interaction Hsa between system and apparatus 
acts only briefiy while entanglement of the system and apparatus states is established, but is 
unimportant during the subsequent evolution of the density matrix that results in the actual 
measurement. He also makes the simplifying assumption that the states which actually 
distinguish between quantities being measured have equal eigenvalues of the non-interaction 
part of the Hamiltonian Hs + Ha + H^, which implies that for the submatrix of spanned 
by these states, the commutator [Hs+Ha + Hs, E[p]\ is zero, and so these commutator terms 
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in Eq. (24b) can be neglected. With these simphfications, Eq. (24b) becomes 

^ = -i[HAE. E\p\\ - \AHm. [Has, E[p]]] , (30a) 

or when the non-Schrodinger term is omitted, as in Zurek's analysis, 

^ = -^[H^,E[p]] . (306) 

Zurek points out that the evolution of Eq. (30b) introduces correlations between the appara- 
tus and the environment, which select as the "pointer basis" of the apparatus, that registers 
the measurement, the eigenstates \Ap) of a "pointer observable" ft that commutes with H^; 
in other words, the pointer basis projectors Up = \Ap){Ap\ must satisfy 

[Up,H^]=0 . (30c) 

Returning to the full evolution equation of Eq. (30a), with the non-Schrodinger terms in- 
cluded, we see that the argument of Sec. Ill, when applied to this equation using Eq. (30c), 
imphes state vector collapse to the eigenstates of the Zurek pointer basis. Thus an energy- 
based stochastic reduction equation, when analyzed within the framework of Zurek's approx- 
imations, is consistent with, and adds further support to, the picture of the measurement 
process that Zurek proposes in [18]. 

In addition to the issues just discussed, there are further questions that must be 
addressed in an energy-based approach, such as whether Hughston's estimated a gives suffi- 
ciently rapid (but also not too rapid) reduction of state vectors for all classes of experiments 
that have been carried out. Answering this question is beyond the scope of the present pa- 
per, but is an important issue for future study. Ultimately, the decision between an energy- 
based or localization-based approach (or yet some other choice of the operator A driving the 
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stochastic terms) may depend on which form of the modified Schrodinger equation can be 
derived as an approximation to relativistically invariant physics at a deeper level. 

To summarize, we have shown that Hughston's stochastic extension of the Schrodinger 
equation has properties that make it a viable physical model for state vector reduction. This 
opens the challenge of seeing whether it can be derived as a phenomenological approximation 
to a fundamental pre-quantum dynamics, along the lines of existing work on open dynami- 
cal systems [19]. Specifically, we suggest that since Adler and Millard [20] have argued that 
quantum mechanics can emerge as the thermodynamics of an underlying non-commutative 
operator dynamics, and since the corrections to the thermodynamic approximation in this 
dynamics are driven by the trace of the energy operator multiplied by a coefficient parame- 
ter with dimensions of inverse mass, it may be possible to show that Hughston's stochastic 
process is the leading statistical fiuctuation correction to this thermodynamics. 
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